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Some exact identities connecting the one- and two-particle Green's functions in the presence of 
spin-orbit coupling have been derived. These identities is similar to the usual Ward identity in the 
particle or charge transport theory and a satisfying spin transport theory in spin-orbit coupling 
system should also preserve these identities. 
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I. INTRODUCTION 

In usual quantum transport theory of charge (or particles), the derivatives of exact single- and two-particle Green's 
functions (or equivalent vertexes) satisfy an exact identity, called generalized Ward identity^ which essentially rep- 
resents the conservation law or gauge invariance. On the other hand, in a practical calculation of the quantum 
transport phenomena we often adopt some approximation methods. Therefore, it is a very important criterion to 
ensure that the approximation remain this exact identity, this kind of approximations is called as self-consistent or 
conserving approximation. 2 It turns out that the self-energy must be taken as a ^-derivative form4 Recently, there 
has been increasing interest in spin transport in spin-orbit coupling systems and some resulting effects, e.g., spin-Hall 
J]3 ■ effect (SHE) ^^^L ^V ! 11 ! 12 ! 13 ! 14 ! 15 ! 16 because it maybe provide a promising method to manipulate the electronic spin 
without application of magnetic fiel d 17 ' 18 ! 19 ' 20 in spintronicsi 21 ' 22 As well known, the spin is not a conserved quantity 
in the presence of spin-orbit coupling, therefore, a natural question is that whether there are also some exact identities 
similar to the usual Ward identity in the charge transport theory? In this Letter we will derive some exact identities 
connecting the one- and two-body Green's functions in the presence of spin-orbit coupling. It is a consequence of a 
' more general continuity-like equation, which was given as an equation valid for average values of the spin-density, 
current and torque in our earlier paper^ here we shall prove that it is valid not only for the averages of these 
quantities, but also holds as an operator identity. Moreover, it will imply some other exact identities connecting 
various correlation functions, which represent more general matrix elements of the spin-density current and torque. 
■ Therefore, a self-consistent or conserving approximation for the response of the spin- density, current and torque must 
£NJ ' also preserve these identities. 
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II. THE MAIN IDENTITIES 
A. Model and Hamiltonian 

Without loss of generality the Hamiltonian with a general spin-orbit coupling can be written as 



^ : fi = / ^(x)Mx,V)] Q/3 V^(x)+ / «(x-x>i(x)^(x')^(x>a(x). (1) 

o ' 
o 



Here and hereafter a summation over repeated suffix is implied and 

h = tt^HV) 2 +e<f> + h so (2) 
Im 

with Vj = di — di — i-A$, i — x,y, z, h so is the spin-orbit coupling term. We suppose 

h so = £W + H^i-id,) + H^i-idiX-idj) + Hgli-idiX-id^-idk) + • ■ • , (3) 

with iJj™?., (n = 0, 1, 2, • • • ) an operator acting on spin variables, meanwhile, it is also the component of a tensor of 
rank n'th and is supposed to be symmetric under the permutations of its suffixes it, • • • ,i n - Moreover, we assume 
that all H^"}. in are independent of x except H^. 
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B. Two Mathematical Lemmas 

Firstly, we have to introduce two useful identities. The first one can be immediately verified by using the technique 
due to Takahashi, Umezawa and Lurie^ 
Lemma. 1 Let 



*(x) = 



, $ t (x) = (0 t 1 (x),4(x),-..^] v ( J 



are functions of x. A^^...^ is a matrix of N x N, and is assumed to be symmetric with respect to the permutations 
of its suffixes. It is easy to verify that for any integer n > 1, we have 



&A ili2 ... in (d? li2 ... in V) - (-l)"(9? li2 ... in $t) Aili2 ... in * = di Ji(x) = V • J(x), 



(4) 



where 



n — 1 n — 1 

= Ec- 1 )*^"^^)^"-*"-^"^-!*) = E(- 1 ) fc$t H---H A M 1 --^- 1 %^.L_ 1 *- ( 5 ) 



fc=0 



fc=0 



Here and hereafter we use 9 to indicate the differentiation acting on its right-hand side factor \f r , while 9 acts on 
its left-hand side factor 

Moreover, we can verify the following theorem: 

Lemma.2 Let n > 1 and 9, = |( 9 j — 9 we have 



where 



^^...^(ar^...^*) - 'i- a, , 97 li2 ...,, * = s 4 j;(x) = v • j'(x), 

(-l) n (flS <a ..^* t )A il i a ... iB * - ^A,,,,..,^,,..^ - 9,Jf (x) = V • J"(x), 



J?W = 5fE $tA -----^t-H^-l- 1 *)=^'( x )» 

\k=o ) 



(6) 



(7) 



and d. 



= 5(9^ - 9 il )---|(9 ira - 9 in ). 



Proof. From the Def.J7|), we can directly verify that 



9, J?(x) - ~ ( 9 4 + 9i) P/(x) = (9i - 9 4 )P,(x) = •!> A , fe ("9£ 



- 9™ ) * 



(8) 



The second line of Eq. ([5]) can be verified in a same way. □ 

The Eq. (j4j and (|6|) also remain valid if all di 's are replaced by 9i 's because 9i = 9 a + i ~ ^ , so 9 j + 9 < = 9 j + 9 , 



C. Generalized equation of continuity 

Now we consider a density operator p a u (x) = (x)o-^g^g(x) = 5' 1 '(x)(t m ^'(x), where 



p a n will represent the particle density if the matrix is a unit matrix er = 
to the Pauli matrix a 1 , er 2 , a 3 . According to the equation of motion 
d 



1 
1 



or the spin density if a 11 equals 



dt 
9 



#t(xi) = 



-* ( x t ) = -i 



* T (x<)/i (x, -9)+ / i;(x-x)# T (xt)# T (xi)#(xt) 
fro(x,~9*)*(xi) + y u(x-x)* t (xf)1'(xi)*(x<) 



(9) 



where the bar over x is to indicate the integral variable, we have 
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j- fPa . (x) = -i*t(xt) ^ (x, 9 ) - ho(x, - 9 Kj *(xt). 
The four-operators terms are canceled because 

y [u(x - x) - u(x - x)]* t (xt)* t (xi)*(xi)*(x<) = 0, 
Substitute Eq.® into Eq.([i"0]) we obtain 

^p CT ,(x) = -i*+(xt) f ^- (aV - <t"V 2 ) + [a^ so (x, 9) - A so (x,-9K \ *(xt) 



9i f 



2 m 



(10) 



(11) 



The right-hand side of the Eq. (Ilip includes two different kinds of terms. The first one comes from the kinetic energy 



2m 



* f (xi) 



A 2 ct ai - cr^V 2 



*(xt) = -V • Jff 



(12) 



where 



J^(xt) = — [*t(xi)o^W(x*) - V^xty^xt)] 



(13) 



The second term is due to the spin-orbit coupling 



*»ho(x, 9) - ft so (x, - 9 )<rH * = -^(-i)"- 1 ff[(7''ffj!!, B 9»... i „ - (-irA ( ;l„^..,J*. (14) 



Using the identities 



(n) 



# (n) - a» = -{a^ H (n) ■ \--\<7> 1 H (n) 



(15) 



Eq. (Trj} becomes 



a^hsoix, d ) - /i so (x, - 9 )a' J 



* = E^- $t ^'^K^U-(-i) B as^) 



2 

r(») 



+ K,^. i j(9».., n + (-ir97 l ..,.j]*. 



(16) 



Let & — ^ , A^...,^ = {H^.i jff 11 } and use the Lemma. 1, the first term of Ea. pB|) can be written as V- J^, where 



(17) 



fc=0 



The second term of Eq. can be rewritten as V • JJ, — T' 7 ' , where the T a>l and J' 7 '' are respectively defined as 



00 



9 - *. 

l\ ■■■l n 



(18) 
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n-l 



(19) 



n>2 



We notice that the T a ^ is given by a polynomial of <9j = l/2(<9j — <9j), therefore, it can not be rewritten as a 
divergence form3 Therefore, we obtain the following identity 



dp a v (xt) 
dt 



V • 



Jf (xt) + 3f Q (xt) + j£ (xt) = (xt) . 



(20) 



The generalized continuity equation given in our previous paper— can be obtained by takeing the average value on 
both sides of Eq. (l2"01) . Evidently, Eq. (|2"D|) includes two different cases: 1). — <r°, in this case p„o is the number 
density p, from the Def. (fT5|) and Dcf. flHj) we have T a — 0, J^ Q = 0, so, 3 a + J° is the current density J and we 
have the following local conservation law of number 



dp(xt) 
dt 



V-(J (xt)+J so (xt)) = 0. 



(21) 



2). cr^ = a 1 , a 2 , a 3 , Eq. (j20|) represents the continuity- like equation of spin. 



D. Exact identities connecting the correlation functions 

Because Eq. (|20l) is an operator identity, it is valid for arbitrary matrix-element of both sides of this equation. 
Therefore, for any operator O(x't'), we have the following identity: 

| (T (xt)O(x't')]) + V • (T[(jf (xt) + j£ (xt) + j£ (xt))0(x't')]) 

= (T[T^(xt)0(x't')]) +( 5(t-t')([^(xt),0(x't')]), (22) 

where the T is the time-ordering operator and O may equal to, e.g., p a ^ or J ^ etc. Similarly, for arbitrary operators 
A and B, we have 

A(r{p^(l)>l(2)B(3)}) + Vx • (T{[jf (1) + j£(l) + j£ (1)]A(2)B(3)}) - (T{T^(1)A(2)B(3)}) 

= 5(h - t 2 )(T{[p a , (uti), A(a;2ti)]B(3)}> + <5(ti - t 3 )(T{A(2)[p CT( .(a: 1 ti), B(z 3 ti)]}). (23) 
For example, if A(2) = tp a (2) and S(3) = ^(3), Eq.jJSSJI becomes 

^■<r{^(i)^ a (2)^j(3)}) + vi-<r{[jf (i) + j s or ;(i)+jr;(i)]^(2)^(3)}> 

= (T{T^(1)^(2)$(3)}> - <5(1 - 2)^(^(2)^(3)}) + 5(1 - 3)(T{^(2)^(3)]})^ 

= ^{1^(1)^(2)^(3)}} - i [5(1 - 2KG(2,3) - 6(1 - 3)G(2, Z)a»] ap . (24) 

These identities represents the exact relations satisfied by the generalized responses functions of the spin- density, 
current and spin torque to some external perturbations. The physical essential of these identities is that the responses 
of the spin-density, current and torque in non-equilibrium should also obey the generalized continuity equation, since 
this equation is an exact operator identity, which would be satisfied in any situation. Therefore, if we calculate the 
spin-density, current and torque by some approximate methods, then the approximation must also preserve these 
exact identities. 



E. Rashba Model 



Because the expressions of the spin-current and torque are dependent on the concrete Hamiltonian, therefore, the 
form of these identities also depend on Hamiltonian of the system. As an example, we consider a 2D system with 
Rashba spin-orbit coupling 



flR = a((T x p y - (TyP x ). 



(25) 
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According to the definition of Def.([l3]). (fT7| and (fT8|) , we have 
(T^(l)Vv(2)^/(3)]> =a^G w ^(l,2;l + ,3), 

<T[jf (1)^(2)^,(3)]) = [(Vi - Vv)G^„(l,2; l',S)] v=1+ , 

(T[J% (1)^(2)^,(3)]) = | V„, « a* V,) G w „ (1, 2; 1+, 3), 

<T[T" M (1)^(2)^(3)]) = |{(K,^](9 X1 -0*,) - -9, ; )V r) G,„ ) ^(l,2;l',3)} 1 , =1+ , (26) 

where G w „'(l,2; l',2') = (T[^(1)^(2)V^ (l')V^ (2')])- Here we only focus our attention on the identities related 
to spin transport, so, take the <r M = a x , a y and a z respectively and substitute Eq. (|2"6|) into ([24]) . we get 

o% v I ^G^'w (1, 2; 1+ 3) + V ~ [(Vi - Vy ) G vvrl i v >(l, 2; l', 3)] 1/=1+ 1 + afi^ G vvnv > (1, 2; 1+, 3) 

- ia<7*,„ - d x ^G^, v , (1, 2; 1', 3)] v=l+ - i [5(1 - 2)^0(2, 3) - 5(1 - 3)G(2, 3)a% v , (27) 

|^-G Wl/ '(l>2; l + ,3) + V • -L [(Vi - V!0G W ^(1,2; l'^)]^ j - a9 xl G^(l,2; l+,3) 

= iaa*, v [{d yi - ^)G Wl/ (1, 2; l', 3)] 1#=1+ - * [5(1 - 2)<r»G(2, 3) - 5(1 - 3)G(2, 3)^, (28) 

{^ Gw, " (1 ' 2; 1+ ' 3) + V ■ 2L [(Vl ~ Vl ') G W-'(1, 2; 1', 3)] 1/=1+ J 
= -ia{[o% v (d Xl -d^+a^idy, -5, ; )]G Wl/ '(l,2;l',3)} i/=i+ 
- i [5(1 - 2)^G(2, 3) - 5(1 - 3)G(2, 3)a\ v , (29) 



If we define L w „, (1, 2; 1', 2') = G W „<(1, 2; 1', 2')-G w /(l, l')G^(2, 2'), and substitute this definition into Eq.(27]) 
and (|2"9"1) , we shall obtain similar identities connecting L and G, in which all G„ V7] i v i (1, 2; 1', 2')'s are replaced by 
L vvrt ' v i(l, 2; 1', 2')'s, while the contribution of the second term G w ' (1, l')G vv > (2, 2') of L will be vanishing owing to 
the general continuity equation. The L is a general response function, which can also be given by 



■^•qvq'v 1 (lj 2, 1 , 2 ) 

where the U dependence of G(U) is given by 

iG w (l,l , ;U) = 



&G VV , (2,2' ';[/) 
5t/^(l',l) 



(7=0 



<T[V>„ (1)^ (l')e _i J £ ^'& s '>*J< 2 >*v &)}) 



For the system with other kind of spin-orbit couplings, e.g., Luttinger model, Dresselhaus spin-orbit coupling or 2D 
cubic Rashba model ) 25 i 26 ' 27 i 28 the Hamiltonian maybe include square or cubic terms in momentum. According to the 
definition, the expressions of the spin-current and torque will become more complex, so, the concrete form of these 
identities will be also become more complex, but the principle is the same. 

III. SUMMARY 

Some exact identities connecting the correlation functions and the one-particle Green's function in spin-orbit cou- 
pling systems were derived. These identities are very similar to the usual Ward identity in charge or particle transport 
theory, although the former are originated from a more general continuity-like equation while the latter is due to the 
conservation law or gauge invariance. These identities can provide a criterion of whether an approximate calcula- 
tion of spin current and spin torque in the presence of spin-orbit coupling is self-consistent. Therefore, in practical 
calculations we must adopt some satisfactory approximations, which can preserve these identities. 
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